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On the Largest Singular Values of Random Matrices with 
Independent Cauchy Entries 

Alexander Soshnikov*and Yan V. Fyodorov''' 



Abstract 



We apply the method of determinants to study the distribution of the largest singular values 
of large m x n real rectangular random matrices with independent Cauchy entries. We show 
that statistical properties of the (rescaled by a factor ^l^^^i ) largest singular values agree in the 
Pm ' limit with the statistics of the Poisson random point process with the intensity —x~^^'^ and, 

I therefore, are different from the Tracy- Widom law. Among other corollaries of our method we 

<~i ' show an interesting connection between the mathematical expectations of the determinants of the 

, complex rectangular mxn standard Wishart ensemble and the real rectangular 2m x 2n standard 

Wishart ensemble. 

m ■ 1 Introduction and Formulation of Results. 
> 

I The main goal of this paper is to study the spectral properties of a large random matrix with i.i.d. 

Tij- ' Cauchy entries. In other words we consider a rectangular mxn matrix A = (uij), 1 < i < m, 1<J< 

. n, where {aij} are independent identically distributed Cauchy random variables with the probability 

^ I density f{x) = ;;^Yq3py- Our goal is to study the singular values of A as the dimensions of a matrix 

O I go to infinity, m ^ oo, n — > cxd. This is clearly equivalent to studying the eigenvalues of a positive- 

c~| ■ definite n x n matrix M = A^A. Matrices of such type are quite often called sample covariance 

, matrices in random matrix literature. Positive-definite matrices are of particular importance in 

^ I statistics (we refer to [^, [22 for the classical works on statistical applications of the spectral 

properties of Wishart matrices, and 22 for a few recent developments and applications to various 
fields). They also are of long-standing interest in nuclear physics, starting with the classical works 
^ I [SO] , ^0] . More recently they were used to model the " dissipative" part of the effective Hamiltonian 

■ in quantum chaotic scattering (see JH] and references therein) and appeared to be very intimately 

connected with the "chiral" ensembles studied in Quantum Chromodynamics, see UHl- As other 
important applications of random positive-definite matrices we mention that they are used in a 
branch of condensed matter theory known as mesoscopics to model famous universal conductance 
fluctuations and other transport properties of small metallic samples and quantum dots, see [1], and 
also emerged in theory of information communication in random environment j45j . 
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It is well known that if the entries of A are i.i.d. random variables with zero mean and finite 
variance o"^, the empirical distribution function of the eigenvalues of ^-A^A. converges in the limit 
m — > oo, n — > oo, m/n ^ 7 G (0,+oo) to the Marchenko-Pastur law (see e.g. [301, PI) defined by 
its density 



= (27rx7CJ^) ^a/(6 — — a), a <x <h, (1) 

where a = ct^(1 — 7""*^^^)^ and b = cr^(l + 7""*^^^)^ (we assume here 7 > 1). Since the spectrum of 
A^A differs from the spectrum of AA^ only by the multiplicity of the eigenvalue A = (for m > n the 
matrix AA^ has m — n additional zero eigenvalues) for the rest of the paper we can assume m > n. 
Under the assumption that the fourth moment of Uij is finite Yin , Bai and Krishnaiah ( , see also 
[35] . 0) showed that the largest eigenvalue of ^A^A converges to b almost surely. Recently Johnstone 
proved that in the standard Wishart case (i.e. {aij} are i.i.d. A^(0, 1) random variables) the properly 
rescaled largest eigenvalue converges in distribution to the /? = 1 Tracy- Widom distribution (( 43 ), 
see also (^]). Soshnikov ([SHI) generalized the result of Johnstone to the non-Gaussian case provided 
n — m = 0(n^/^) and the moments of the matrix entries do not grow very fast. There are quite a few 
standard methods that have been successfully used for Wigner and sample covariance matrices in the 
case when second and higher moments of matrix entries exist, most notably the method of moments 
(UHl, HO], ISO], ISZI, ISHl, inni), the method of resolvents [20], EI), the method of orthogonal 

polynomials (d], HH])) the method used by Johansson (and recently extended by Ben Arous and 
Peche to the case of sample covariance matrices) which is based on the Kazakov-Brezin-Hikami trick 
(l28|,IH|) inD; etc. Unfortunately, the above mentioned approaches are not suitable for the Cauchy 
case. In particular, one can expect the spectral properties of A^A in the Cauchy case to be rather 
different from the case of a finite variance. In our view this makes the studies of the Cauchy case 
especially interesting. Denote the eigenvalues of A^A by Ai > A2 > A3 > . . . > A„. It is expected 
that the majority of the eigenvalues are proportional to mn. We would like to specifically single out 
the reference |12j where Cizeau and Bouchard studied the spectral properties of the Wigner random 
matrices with the heavy tails (see also references |23j for physical papers on the so-called Levy- 
Smirnov unitary ensembles). Among other things Cizeau and Bouchard argued (on a physical level 
of rigor) that the empirical distribution function of the eigenvalues of a properly normalized Wigner 
matrix (with the heavy tails of the marginal distribution of matrix entries) converges to a limiting 
distribution that can be obtained as a solution of a quite complicated system of two integral equations 
(we refer to the formulas (15), (12a), (12b) in |12j). It is not difficult to guess the right order of 
the normalization: after the normalization, the norm of any given matrix raw has to be of order of 
constant (in particular in the Cauchy case one has to normalize the matrix entries by n"^ and in 
the case of a finite variance the normalization is -y/n). The support of the limiting distribution is 
the whole real line. One can expect (arguing at the same level of physical rigor) to derive a similar 
system of integral equations for the limiting distribution function of the eigenvalues of :^A^A in the 
case of the i.i.d. Cauchy entries of A. The support of the limiting distribution should be the positive 
half of the real line. Our results formulated below indicated that the asymptotics of the spectral 
density at infinity should be ^^3/2 ■ This suggests that the largest eigenvalues of A^A grow faster than 
mn, in fact we will show below that the largest eigenvalues are of the order of vn}r? . Let us rescale 
the eigenvalues by that factor: 

= « = l,---,n. (2) 

The goal of this paper is to study the local distribution of the largest eigenvalues by the method of 
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determinants. Our main results are Theorems 1.1 and 1.2 formulated below. 



Theorem 1.1 Let A be a random rectangular m x n matrix (m > n ) with i.i.d. Cauchy entries 
and z a complex number with a positive real part. Then as n ^ oo we have 

— 1/2 ^ / 9 \ 

lim E(det{l + ^-^A'A)) = eY\{1 + zXi)-^/^ = exp ( ^) (3) 

1=1 ^ ^ 

oo 

= ¥.J{{l + zxi)-^/\ (4) 



where we consider the branch of ^fz on D = {z : 5Rz > 0} such that \/T = 1, E denotes the 
mathematical expectation with respect to the random matrix ensemble defined above, E denotes the 
mathematical expectation with respect to the inhomogeneous Poisson random point process on the 
positive half-axis with the intensity ■, o,i^d the convergence is uniform inside D (i.e. it is unform 
on the compact subsets of D). For a real positive z = t^, t G M, one can estimate the rate of 
convergence, namely 

/ f2 \ -1/2 n . s 

lim ^ det(l + ^^>l*^) =£;TT(l+t2A0-i/2 = exp — |f|(l + o(n-i/2+^)) , (5) 

n-^oo V m'^n'^ J V vr / 

where e is an arbitrary small positive number and the convergence is uniform on the compact subsets 
o/[0,+cx)). 

We will discuss the properties of Poisson random point processes at the Appendix. A very useful 
introduction to the elementary theory and methods of random point processes is ^Sj. It is not a 
coincidence that the intensity of the Poisson random process in the above theorem is equal to the 
leading term of the asymptotics of the density of the square of a standard Cauchy random variable. 

We claim that the result can be generalized to the case of a sparse random matrix with Cauchy 
entries. Let as before {a^fc}, 1 < i < m, 1 < /c < n, be i.i.d. Cauchy random variables, and 
Q = {ijk) be a m X n non-random rectangular — 1 matrix such that the number of non-zero entries 
in each column is fixed and equals to 6„. For technical reasons we assume that b^ grows to infinity 
as some power of 71, i.e. bji ^ 71^, for some ^ ^ 1 and ln?Ti is much smaller than any power of 
n. We define a m x n rectangular matrix A with the entries Tjk = QjkO-jk^ ^ ^ j ^ inn, 1 < k < n. 
As before we denote by Ai > A2 . . . > A„ the eigenvalues of r*r. The appropriate rescaling for the 
largest eigenvalues in this case is going to be Aj = ^^'^^ , i = 1, . . . ,n. We claim that the result of 
the Theorem 1.1 can be extended to the case of a sparse random matrix T. 

Theorem 1.2 Let T be a sparse random rectangular m x n matrix (m> n ) defined as above and 
z a complex number with a positive real part. Then as n —> 00 we have 

lim.E;(det(l + --|^r*r)J = e\1{1 + z\)-^'^ =exY>[--^j (6) 



n^oo y ■ rn^b"^ 



E^(l + zx,)-l/^ (7) 



i=l 
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where, as in Theorem 1.1, we consider the branch of \[z on L> = {z : > 0} such that Vl = I, E 
denotes the mathematical expectation with respect to the random matrix ensemble defined in the para- 
graph above the theorem, E denotes the mathematical expectation with respect to the inhomogeneous 
Poisson random point process on the positive half-axis with the intensity 1 ^^i^d the convergence 
is uniform inside D (i.e. it is unform on the compact subsets of D). For a real positive z = t^, t G M, 
one can get an estimate on the rate of convergence, namely 

f f2 \ "1/2 n , . 

E (^det(i + ;^r*r)J = ^n^^ + = '''^^ ^ o(^n'/'+^)) J , (8) 

where e is an arbitrary small positive number and the convergence is uniform on the compact subsets 
o/[0,+oo). 

The result of Theorem 1.2 can be generahzed even further. Let the setting be as in Theorem 1.2 but 
relax the condition that the number of non-zero entries in each column is exactly 6„ to the condition 
X]fc=i Qjk = + o(l)), j = 1, . . . ,m, (for example the relaxed condition is satisfied by a typical 
realization of a random matrix Q with independent Bernoulli 0—1 entries with Pr{qjk = 1) = bn/n). 

Then we still have E ( det(l + ^^^|p-r*r) j = exp (— -y^) (l + o(l)). The proof is almost identical 
to the proof of Theorem 1.2 and will be left to the reader. 

The case when the number of non-zero terms of Q is fixed in each raw can be treated in a similar 
manner. 

An important consequence of Theorems 1.1 and 1.2 is that the statistical properties of the largest 
eigenvalues Ai, A2, . . . are drastically different from the statistical properties of the (rescaled) largest 
eigenvalues in GOE ([ISl) and real Wishart case (|26i) that are described by the (/3 = 1) Tracy- Widom 
law. 

Theorems 1.1 and 1.2 follow from the Proposition 1 formulated below. 

Proposition 1 Let A = (aji.) be a random rectangular m x n matrix with independent (not neces- 
sarily identically distributed) entries with the characteristic functions of the matrix entries gjk{s) = 
E exp(isajii.). Let ij > 0, i = 1, . . . ,r be some positive parameters. Then the following formula holds 

/ r \ -1/2 

E Yl det(l + tfA'A) = (27r)-'^("+"^)/2 ^ 

/ nd".«d>« exp('-^(|s»|V2 + |p«lV2)) n 9,k(j2'^pfs^^](9) 

Theorems 1.1 and 1.2 imply several important corollaries given in the next section. 

The proofs of Theorems 1.1 and 1.2 and Proposition 1 are given in the next section. Section 3 is 
devoted to application of the method of determinants to random matrices with i.i.d. complex entries. 
We prove in section 3 the analogue of Proposition 1 (Proposition 2) in the complex case and as a 
corollary establish an interesting connection between the determinants in the 2m x 2n rectangular 
real Wishart case and m x n rectangular complex Wishart case (see Lemma 1). 
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2 Proofs of Theorems 1.1 and 1.2 

We start with the proof of the Proposition 1. Consider (det(l + t'^A^A)) , t > 0. Let s 
(si, . . . , s„)*, p = {pi, . . . ,pmY be real n— and m— dimensional column vectors. Let B{t) 

Then 

-1/2 



(n+m)/2 



J d^s d'^p exp{-{s,p)B{t)is,pY) 

j d^'s d!^p exp(-(|s|2 + exp(-2i ^ ta.fcp^Sfc) j . (10) 

\ l<i<m,l<fc<n / 



^ ^ {n+m)/2 

vr . 



In pop we used standard properties of the Gaussian integral. The formula (jlOl) and the indepen- 
dence of matrix entries imply 



fr \-i/2 
S mdet(l + t2^*.4) 



vr 



\j=i 

-r(7i+rrt)/2 



/ nd«.«d>« exp('-^(i.«iv2+ip«iv2)) n 5. 4^^)11) 

The Proposition is proven. 

To prove Theorem 1.1 we observe that the functions at the l.h.s. of Q are analytic and uniformly 
bounded in D = {z : 5Rz > 0}. Therefore, by the Vitali's theorem it is enough to prove the 
convergence for real positive z. Let us denote z = t^^ where t is a positive real number and apply 
the result of the proposition in the case r = 1. Since the matrix entries of A are i.i.d. Cauchy we 
have gjk{s) = g{s) = exp(— and 

E{det{l + t^A'A))~^^^ = 

(2^)-{n+m)/2 /■ exp(-(|sp + |p|2)/2) X TT exp (-t | SfcPj | ) = 

■^-^ l<j<n,l<A;<m 

2- / d".(27r)-"/2exp('-l|.|2') / (2vr)-"/2 exp [ -i f;(p| + 2p,t ) 

^ ^ •^^+ \ j=i k=i J 
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= / d^si27r)-^/''exp(-l-\s\A^"'{tj2\sk\), (12) 

where = 2e2^'/2 j+oo _^g-_tV2^^^ particular, ^'(0) = 1 and ^-'(0) = -(2/7r)V2. it is easy 

to see that the function ^{y) is monotonly decreasing on [0, +oo), in particular 1 = ^'(0) = 
max[o^+oo) ^(y)- Indeed, ^'(y) = 2ye^ ^"^ Iy°° ''^^^ ~ '^^^ assertion then follows from 

the inequality y~^e~^^/^ > /^^°° e'^^/'^dt for y > 0. 
Replacing t by ^ we arrive at 

The r.h.s. of the last formula suggests to use the law of large numbers and large deviations estimates 
for the sum of the absolute values of n standard Gaussian random variables. Since 
(2/7r)V2, we see that lYTj=i\sj\ = (2/vr)^/2 + ^(^-1/2+^) ^j^j^ probability 1 - 0(exp(- 



Recalling that ^'(0) = 1 and ^-'(0) = -(2/7r)i/2 we get 

S(de.(l + J%.4'A,)"'^ / d",(2,)-"/=expf-i|»p) / _ 2.(1 + o(„-/^^.)) V" 
\ n'^m'^ J Jjin V 2 / \ vrm / 



2,(1+ o(n-V^^-)r^ / 2.(1 + o(,,-/'t-)\ 
\ 7rm / \ TT / 

for any e > 0. Theorem 1.1 is then follows by the Vitali theorem. 

The proof of Theorem 1.2 is very similar. Again we can restrict our attention to the case when 
z is a real positive number, z = t^. We have 

£;(det(i + t2r*r))"^/^ = 

\ m n 

-^isnU'^itY.qjkiski). (15) 
^ j=i k=i 

Let k[^\ : ■ ■ ■ , be the indices k for which qjk = 1. Then X^^^i qjk\sk\ = Yl'i=i \^i.U) \ and we 



can claim that ^ X]f=i I^.O)! = (2/vr)^/^ + o(6„^''^+') with probability 1 - 0{exp{-bl;^^^) for each 
1 < j < m. Since we assumed that ln(m) is much smaller than any power of n we get similarly to 
(dH) that 



m 



fl-^'"+°'^"""''y" = expf-^'"+''<^"""''V (16) 
Theorem 1.2 is proven. 
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Below we restrict our attention to the corollaries of Theorem 1.1 (full matrix case). In corollaries 
of Theorem 1.2 are basically identical to those of Theorem 1.1 (with an obvious change of n to bn 
where it is needed). 

Remark We are not aware that the results of Theorem 1.1 and 1.2 are enough to imply that 
ths statistics of the largest eigenvalues are Poisson in the limit of n oo. Indeed, to prove the 
Poisson statistics in the limit one has to show that 

n +00 

\mi^El[{l + f{X,))=Ell{l + f{x^)) (17) 

^ i=l i=l 

for a sufficiently large class of the test functions /. The results of Theorems 1.1-1.2 claim that (|17|) 
is valid for f{x) = (1 + zx)~^^'^ — 1 for all z such that > 0. Below we formulate several corollaries 
of our main result that are weaker than the claim about the Poisson statistics, but still give us some 
information about the behavior of the largest eigenvalues. The proof of the Poisson statistics for the 
largest eigenvalues in Wigner random matrices with heavy tails will appear in |4nj . 
Remark 1 

It follows immediately from the result of the Theorem 1.1 that "only a finite number" of the 
eigenvalues Aj are of the order of n^m?. Indeed, let Nn^m be an integer growing to infinity arbitrary 

slow as n ^ CO and > be an arbitraty small positive number. Then Pr^#(Aj > 5n?m?) > 
Nn,rn^ ^ as fi ^ OO. Indeed, suppose this is not the case. Then Pr(^#{Xi > 8n^m^) > iV„,„) > 
K > and E (det(l + ) < k(1 + t^K)"^"-" + (1 - k)1. One obtains a contradiction 



since for t = N^l^ the r.h.s. of the last inequality does not go to zero (see also next remark). One 
can also rewrite the statement of this remark in the following way: for any positive 5, k there exist 
no(5, k) and C{5, k) such that Pr ^#(Ai > dn^rn^) > < k for all uq < n < m. 
Remark 2 

It is clear from the proof of the Theorem 1.1 that the asymptotic result 

t2 I 2t(l+o(n-l/2+e))\ 



/ ^2 

E det(l + -^^A^A) = exp 



vr 



holds uniformly in t on compact subsets of [0,+oo). In particular the result is valid for a sequence 

Below we formulate and prove some additional consequences of Theorem 1.1. Our first observation 
is that one can repeatedly differentiate (jHI with respect to parameter z. 

Corollary 1 iei Sftz > 0. Then 

hm eT\{\^z\,)-^I^ I I =lz-V2expf--^') = (18) 



CO 



where, as above, E stands for the mathematical expectation with respect to the inhomogeneous Poisson 
random point process on (0, +oo) with the intensity p{x) = ^^3/2 ■ 
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Remark 3 

If we let z ^ in (|18|) one gets E (^^^i A,^ = +00, which trivially follows from the fact that 
matrix entries of A are Cauchy random variables. Essentially the result of the corollary can tell us 
how fast the mathematical expectation E ( Xj j grows if we restrict the summation only to Aj < L 



where L is large. 

Proof of the Corollary 1. The result immediately follows form the uniform convergence of the 
analytic functions in 1)314(1 . 

By differentiating (jSJ twice one can obtain in a similar fashion that 



lim EY\{l + zXi) 

i=l 



4 _i 2 



.^-1/2 



vr^ vr 




Ji^-i + 1^-3/2 ] I --^ ] . (20) 



Corollary 2 There is a contant C which depends on 7 such that for Pr{j^^ > x) < Cx 
uniformly for large n <m and x. 

Indeed, it follows from Theorem 1.1 and Remark 1 that E{1 + t^Ai)~^/^ > exp ( — 2t{l + 
o(n^^/^+'^))/7r) uniformly in t on compact subsets of [0,oo). Therefore (1 — Pr(Ai > x)) + ^-^^^^2^ ^ 
Pr{\i > j;) > exp - |t(l + o(n-i/2+e))^ ^ ^j^jpi^ implies 1 - exp - |t(l + o{n-^/^+^))^ > 

(1 - ^7=5^) Pr{Xi > x). Choosing t^x = 1 we obtain that |7x"^/^(l + 5) > {1 - 2-^/^)Pr{Xi > x) 
for all sufficiently large n, m, and x. 
Remark 4 

It is not difficult to show that in probability Ai = 0{n'^m?). To see this we observe that the 
operator norm ||A|| can be bounded from below by maxi<j<m,i<A,<n The maximum of n x 

m i.i.d. Cauchy random variables is of the order 0{nm) (with the limiting distribution of ^ x 
maxi<j<m,i<fc<n lojfcl easily computable, namely Pr maxi<j<m^i<fc<„ |ajfc| < x) ^ exp(-^)). 
We expect that the limiting distribution ^^^^ also exists and is given by the distribution of the right- 
most particle in the Poisson process with the intensity p{x) = ^x~^/'^ ^ in other words, lim„_+oo Pr(Ai < 
x) = exp(— 

It is a useful excercise to see what Proposition 1 gives in the Wishart case. Below we treat the 
case of one determinant. 

E{dei{l + t^A'A)y^'^ = 

(2^)-{n+-)/2 I d^.^p exp (-(|s|2 + |p|2)/2) X n exp f-htSkPjf) = 

l<j<n,l<k<m V ^ J 

/^^d".(2.)-/^exp(-lN^) /^^0(2.)-W2nexp(^4pi(l + .^|;.,^)j ^ 
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r+oo 1 
= Cn exp(--r2)r"-^(l+iV)-™/2rfr 
^0 2 

[•+00 

= c„2t-M e-Vt-i(l + 2t2r)-'"/2dr, (21) 
Jo 

where c„ is the normahzation constant, = exp(— ir^)r"~^dr = 22~'^T{n/2). To study the 
global distribution of the eigenvalues in the Wishart ensemble one has to consider rescaling ^ ^ 
(since typical eigenvalues of A^A are of the order of n). It follows that 

/ ^2 \ -1/2 /.+00 

S det(l + -^*^) =[r(n/2)]-M e-Vt-i(l + 2t2-)-W2^r = 

V ^ / Jo 

r+oo 

[r(n/2)]-i (n/2)"/2 / g^t Vt (1 + iVj-^/V-^dr = 
Jo 

r+oo 

[r(n/2)]-^ (n/2)"/2 / e'^^^'^-^dz, (22) 



where i2(2;) = -2 + ^ ln(l + t'^z) — Inz. The asymptotics of the last integral can be obtained by the 
steepest descent method. The formulas are especially simple in the square case m = n. One then 
can find a positive solution of the equation 

^ = 1 + -—^ --= ^ , / = (23) 
dz l + t^z z z{l + t^z) ^ ^ 



to be z{t) = "^^^2*^ ~- Taking into account that = \/4t^~-i-T at z = z{t) we obtain that in 
the square Wishart case 



t 



2 \ -1/2 



E det(l + -A A) 
n 



2"n"/2 / -1 + V4t2 + l\ ^ ^- ^ „ / 27rV4t2 + 1 2*^ ^ ^ ^ 

_e.p (^-n j (1 + V^)-^ „ _i^Vtf.,i <l + 

The fact that the asymptotics in 1)241) is exponential in n is standard . Indeed, it is a straightforward 
excercise to verify that 

lim i In f £; f det(l + —A*A)\ \ = -- [ (ln(l + t'^x)pi(x)dx, (25) 
n^oo nyv ™ / / ^ Jo 

where pi{x) is the probability density of the Marchenko-Pastur law defined in which reflects the 
law of large numbers for the linear statistics Yl7=i ln(l+t'^Ain~^) where {AjjJL]^ are the eigenvalues of 
the real Wishart matrix. The variance of the linear statistics is bounded and has a limit as n — > oo 
(for a rather general class of polynomial ensembles of random matrices it was first discovered by 
Johansson in ([IHl)) thus contributing a constant term to the r.h.s. of H24|) . so that 

^ (^det(l + ^A*A)^ ^ =E(^exp(^-^J2ln{l + t^Xin-^)Jj 
= exp(^-^E{pM^ + t'X,n-')))+^d{t) + o{l)y (26) 
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where d{t) := lim„^oo ^ar ( ^"^-^(ln(l + t^Ajn"^))) and can be explicitely calculated (see (ESI), 
Theorem 2.4). 
Remark 5 

Another class of random matrices we are particularly interested in are Rademacher random 
matrices (i.e. square random matrices with ±1 i.i.d. entries), which we are going to denote by TZ. 
It appears that the questions of the invertibility of a Rademacher random matrix and the estimate 
of the norm of the inverse are of big importance in geometric functional analysis (for example in 
connection with a deterministic construction of Euclidean sections of convex bodies). Similar to the 
previous analysis one can obtain 

1 1/"" /l" \" 

E — , = - — — / 1 iduidviexT) — > (uf + vf) \ \ cos itUjVk)- (27) 

,/M{TTlmn) (27r)«^2„ii ^[ 2{rC^ ) ,^ti ^ ^ 

The fact that a Rademacher matrix R is invertible with probability going to 1 as n ^ oo was 
proved by J.Komlos (see e.g. [HI, chapter 14). More recently, J.Kahn, J.Komlos and E.Szemeredi 
(123) proved that the probability that R is invertible is exponentially close to 1. To the best of our 
knowledge there is no known estimate on the norm of the inverse matrix (which, in our language, 
corresponds to the estimate of the smallest eigenvalue of R^R). 



3 Complex Matrices with i.i.d. Entries 

In this section we consider the ensemble of n x n complex random matrices M = A* A, A = 
iAjk)i<j<m, i<k<n with the joint distribution of the matrix entries of A given by the formula 

Pv{A)dAdA= Yl {d^ajkd'^ajk-fi\ajk\^)}. (28) 

In other words {ajk, 1 < J < ^n-, 1 < A; < n} are independent indentically distributed random 
variables with a distribution depending only on the radial component, and /(x) is the density of 
the distribution of |ajfcP- The ensemble (|28() is a generalization of the standard Wishart (Laguerre) 
ensemble which corresponds to the choice f{x) = e~^. 

In the standard Wishart (Laguerre) case it is known that the smallest eigenvalues are proportional 
to and the (rescaled) A:-point correlation functions are given in the limit n ^ oo by the 

determinants 

Pk{xi, . . . ,Xk) = det{K{xi,Xj))ij=i^,„^k, /c = 1,2,3, ... , (29) 
with the Bessel kernel (with a = 0). 

^H^^ ^ Ja{2V^)VyJ'a{2^) - U2^)^J'^{2^) ^^^^ 

x-y 

where J^, is the J-Bessel function, appears as the limit of the rescaled correlation kernel at the hard 
edge in the Laguerre and Jacobi ensembles (see e.g. ^H], EH)- 

Ban Arous and Peche (0) following the approach suggested by Johansson i [24] ) for Wigner 
matrices have recently shown universality of the limiting distribution of the smallest eigenvalues (as 
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well as in the bulk of the spectrum) for a special class of sample covariance matrices. Their technique 
requires that entries of A have a Gaussian component. 

One of the possible ways to attack this problem for an ensemble ()28() (assuming that all moments 
exist, i.e. and do not grow very fast) is to study the mathematical expectation of the ratio of 
determinants 

^ det(l + r?MM) 

for appropriately scaled (large) real numbers rji, . . . ,rjk,ti, . . . ,ti. For the standard complex Wishart 
(Laguerre) case these expectation values were calculated exactly for any k,l,n, see ^7], |U and ref- 
erences therein, and also used to address objects interesting in mesoscopic physics ^H]) and Quantum 
Chromodynamics HI . 

In particular, one can easily see that ■^Z{r],t)\ri=t = -E' l^"=i ^/!^\ ; where < Ai < A2 < • • • < 
A„ are the eigenvalues of A* A. Such an object can be used to extract the mean eigenvalue density. 
In a similar fashion, by taking partial derivatives of Z{r]i, . . . , ti, ■ ■ ■ ,ti) of higher orders, one can 
study the correlations of the eigenvalues of A* A. To show the universality of the distribution of the 
smallest eigenvalues we need to show that local statistical quantities at the edge of the spectrum 
(near the origin) do not depend (in the limit n — > 00) on the second and higher moments of /. 

The next proposition is analogous to Proposition 1 in the real case. 

Proposition 2 Let A be a random rectangular m x n matrix with the probability distribution given 
by \2t^) . Let ti > 0, I = 1, . . . ,r be some positive parameters. Then the following formula holds 

E (f\det{l + t'^iA'A)\ =vr-"("+'")/ [ ff d"3f?s(')d"9s(^)(i"^3f?/)d'"9/) 

expf-^(|.(')|^ + |/)|^)) n G(\±t4^^A, (31) 

\ 1=1 / l<j<n,l<k<m \ 1=1 / 

where s*-'^ = (s^^\ . . . , s^n) are n-dimensional complex vectors, p^'-* = (pf\ . . . ,Pm) are m- dimensional 
complex vectors, / = 1, . . . , r, 

1 /■2-iT r+co . , r+00 

G'(?/) = — d9j^ f{x)dxe^p(2i{xy)'/\os{9)) = dxf{x)4>{xy), (32) 

and 

Hx) = ^ij^x' = M2x'/'). (33) 
1=0 ^'-^ 

In the special case of single determinant (r = 1) the formula Hell]) can be simplified 

. n m 

E{det{l + t'^A*A)) = / TT e-^'^dui TT e-^'dv^ TT G {t^um) . (34) 

J{o,+<^r+- i=i j=i k,i 

Remark 6 As in the section 2 we can consider the case when the the matrix entries {|ajfc|} are 
independent but not identically distributed with the densities fjk{\x\), 1 < j < m, 1 < n < n. The 
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result of Proposition 2 still holds true provided we replace G in (|^T|) and H32() by Gjk{y) = ^ x 
So'' d9 Jq°° fjkix)dxexp (2i(xy)^/^ cos(6l)) = f^°^ dxfjk{x)<f>{xy). 
Remark 7 

In Jo{z) = 2^ exp {iz cos 9) d9 is the standard Bessel function (@7j). 
Remark 8 

In the special Wishart (Laguerre) case ( which corresponds to f{x) = e~^) one has G{y) = e~^. 
If all moments of f{x) exist and do not grow very fast one can write G{y) = ^^t^ ^"pjr^y'j where 
{on.}n>i are the moments of f{x). 

Proof of Proposition 2. 

Let s = (si, . . . , Sn)*, p = {pi, . . . ,PmY be complex n— and m— dimensional column vectors 
and s* = (sT, . . . , s^), p* = {pT, ■ ■ ■ iPm)- In what follows d^^s and d'^"^p will stand for d'^^sd^Qs 
and d^^pd^Qp correspondingly. Then 

(det(l + t^^M))-^=(det(^^^. *f ))" 

(1 \ n+m 

-j J d^^'s d^'^p exp ( - {s\p*)B{t){s,pf) 

-j / d^"s d^^'p exp(-(|sp + |pp))exp -z ^ (tajkPkSj + tajj^s 
' \ i<i<",i<fc<m 

where as before = ( J and d^ng^^am^ ^ j-jn^^ ^syj^.^cj^. syj^.j^^^ 

We can then write down 

E(f[detil + tjA*A)\ = 



TT 



\Z=1 

-r(n+m) 



/" fTd2"s«d2n^^« exp f - ^ (| s« 1 2 + |/) |2)^ 



n exp -iiauj tisfpf + ar^Y.tisfpf) (35) 

l<i<™,l<fc<m \ \ /=1 (=1 / / 



Let u = tis'J^p^[' . Then we can write 



(exp(— z(ajfcM + ajyfcu)) = — / (i5Rzd9z/(|zp) exp(— i(zu + zn)) 



TT 

f27r /'+00 



1 /"Ztt r+oo 

- d9 drrf{r'^)exp(-ir{e'% + e-'^u)) 
^ Jo Jo 

- / / dx/(x)exp(-iV^(e'^^i + e-*^u)) 
!^ Jo Jo 



27r 

I 

^0 

dxf{x)(f){x\u\'"^), (36) 



27r 

+00 

2 
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where (pix) has been defined in H33|) . Combining (|35() and H36() we arrive at 

^n(det(l + tM*yl))-' = / M / ^2n^(/) rf2,n^(0 gxp (- ^ | | ^ - ^ |2)) 

1=1 «=1'^ 1=1 fc=l 



n GiiY^us^pp^^n. (37) 

l<i<n,l<fe<'Tt V 1=1 / 

In the special case r = 1 the formula can be simplified further 

(1 \ m+n n m 

- ] / cP'^s (f^p exp(- V |s,f - V X 

TtJ jR2{n + m) .^^ 

l<i<",l<fc<m ° I 1=0 ^ ) 

,. n m In m \ n m 

= 2-+" / n ^^^^^ n p^dp, exp - - E n n ^ i^'^ip") 

Ao,+-r+"f=i ,ti \^ ^ ) t=iT=i 

I, n m 

= / ne-"»dn,TTe-'^»dt;iTTG(tVt'/). (38) 

ho,+^)^+- t=i f=i tt 

In the Wishart case one can simplify things even further. The calculations are very similar to 
the real Wishart case considered in Remark 4 (section 2). Since typical eigenvalues of A* A are 
of the order of n the scaling of by a factor 1/n allows us to study the limiting distribution 
of the eigenvalues. It follows from ()38() that we are left with the task of evaluating the integral 
£;(det(l + '^A*A)y' nr=i e-'^'dui Y[%, e^^dv, exp (-f ELi YZi v) , which can be 

reduced to 



L 



e "^'dui e "'dvj exp ^ ^ z;; 

i=l j=l \ fc=l «=1 / 

/ 2 " \-"^ 



(0,+oo)" .^^ 

7 TTT^" 1 + 1 dz 

[n — 1)1 V 

„7i /■+00 /■+00 



r+co n r+oa 

, / ^"-^e— (1 + 1^^)-'" dz = ^ e-^^^^h-'dz (39) 

r(n) 7o r(n) 7o 

where £(2:) = 2: + ^ ln(l + t'^z) — Inz. It is remarkable that the formulas in the complex case are 
identical to those in the real case ()2H22|) modulo trivial change of parameters. We thus proved the 
following result. 



Lemma 1 Let m and n he positive integers and z G C \ (— cxd,0). Then 

z 
2 



E2m. ,2n,real 

(det(l + ^A'A)) = Em ,n,complex 

{detil + zA*A))-^ (40) 
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where at the l.h.s. we have the mathematical expectation with respect to the ensemble of rectangular 
2m X 2n real matrices A with i.i.d. standard Gaussian entries (standard real Wishart ensemble), and 
the r.h.s. we have the mathematical expectation with respect to the ensemble of rectangular m x n 
complex matrices A with i.i.d. standard Gaussian entries (standard complex Wishart ensemble) 

As before, it was enough to prove the result for positive real z = t'^. We remind the reader that 
in the standard real Wishart case all entries {oj^fc} are i.i.d. N(0,1) random variables, and in the 
standard complex Wishart case all entries {SRaj^fc; 9oj,fc} are i.i.d. N(0, 1/2) random variables (so in 
both cases E\aj^k\'^ = !)• 

Remark 9 

If all moments of f{x) exist then under some technical conditions the asymptotics of E (det (l + t'^A*A^^ 
in the global regime depends on the first moment and second moments of f{x)dx, (i.e. on the second 
and fourth moments of the matrix entries Aj^i. This phenomena is known in random matrix the- 
ory: for example in the case of Wigner random matrices, the limiting distribution of a global linear 
statistics Tr h(A) — E[Tr h{A)) (where /i is a test function, say a polynomial, and A is a random 
Wigner matrix normalized so that a typical eigenvalue is of the order of a constant) depends on the 
second and fourth moments of the matrix entries (see e.g. [211 5 IH^)- It is conjectured (and in a few 
interesting special cases verified) that in the local regime the dependence on the fourth moment goes 
away. 



4 Appendix 

A Poisson random point process on the positive half-axis with the locally integrable intensity func- 
tion p{x) is defined in such a way that the counting functions (e.g. numbers of particles) in the 
disjoint intervals /i, . . . , are independent Poison random variables with the parameters Jj p{x) x 
dx, j = 1, . . . ,k. Equavalently, one can define the Poisson random point by requiring that the /c-point 
correlations functions are given by the products of one-point correlation functions (intensities), i.e. 
Pk{xi,...,Xk) = U'j=iPixj). 

Let / : (0, +00) — > C be a test function with a nice behavior at the origin and infinity. Then 

00 00 k OD ^ „ k 

Ej];(l + /(xi)) = l + ^E Yl n^(^^.)=Efc!/ ^l{f{xj)pk{xu...,Xk)dxi---dxk 



X] 77 ( / f{x)p{x)dx ) = exp ( / f{x)p{x)dx ) (41) 
^' ^ '(0,+oo) / \J(o,+oo) J 



k\ 

k=0 



If the test function /(x) equals (1 + zx) ^1"^ — 1, and p{x) = ^^^i^ we have 



1 2 

f{x)p(x)dx= I {{1 + zx)-^/^ - l)—-dx = 



(0,+oo) J{0,+oo) 



which is exactly the exponent in (jlj. For random Schrodinger operators the Poisson statistics of the 
eigenvalues in the localization regime was first proved by Molchanov in [33j (see also j32] ) . 

Acknowledgements. It is a pleasure to thank Boris Khoruzhenko for very useful comments. 
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